Suppression of decoherence effects in the quantum kicked rotor 
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We describe a method allowing transient suppression of decoherence effects on the atom-optics 
realization of the kicked rotor. The system is prepared in an initial state with a momentum distribu- 
tion concentrated in an interval much sharper than the Brillouin zone; the measure of the momentum 
distribution is restricted to this interval of quasimomenta: As most of the atoms undergoing deco- 
herence processes fall outside this detection range and thus are not detected, the measured signal is 
effectively decoherence-free. 
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The kicked rotor has been studied for many years in 
the helm of both classical and quantum Hamiltonian dy- 
namics, where it has the status of a paradigm. Its for- 
mal simplicity, and the richness of its dynamics make it 
one of the most studied dynamical systems. Experiments 
started with the observation by Raizen and co-workers 
[H of quantum chaos in an atom-optics realization of the 
kicked rotor. An impressive number of works are devoted 
to different aspects of the kicked rotor dynamics, includ- 
ing dynamical localization [1 6j , quantum resonances 0- 
llll . sub-Fourier resonances |l2 , chaos-assisted tunneling 
|l3l- 15|. ratchets [l6l - [ll |. and, recently, the observation of 
the Anderson metal-insulator transition [6, 19, 2(|]. The 
atomic kicked rotor is a very clean and flexible system, 
but it is limited by the unavoidable presence of deco- 
herence, mainly due to spontaneous emission 0, Si. In 
the present work, we propose and analyze a method for 
suppressing to a large amount, if only transiently, deco- 
herence effects in the kicked rotor's dynamics, which is 
directly applicable to existing experimental setups. 

Consider a two-level atom placed in a standing wave 
generated by two counterpropagating beams of intensity 
/ and wavenumber k^. On the one hand, the atom may 
interact non-dissipatively with the radiation by absorb- 
ing a photon in one of the beams and emitting a pho- 
ton in the other beam by stimulated emission. This is 
a conservative elementary process in which a "quantum" 
of momentum 2hki^ is exchanged between the atom and 
the field. It generates a potential acting on the center 
of mass of the atom, which turns out to be proportional 
to r2(a;)^/Ai~ (//Al) cos(2A:lx), where fl is the (res- 
onant) Rabi frequency, is the laser-atom detuning, 
(Ai ^ r where F is the natural width of the excited 
state) [H, 22|. This "optical potential" only couples mo- 
mentum states separated by multiples of therefore, 
quasimomentum is conserved. On the other hand, there 
can be spontaneous emission (SE) process which are dis- 
sipative and thus are fatal to quantum effects. The total 
momentum exchange (along the direction of the standing 
wave) in a process involving SE is hkL{l + cos9) where 
6* is a random angle in the range [0,27r). After such a 
process the atom momentum distribution is not anymore 



concentrated around multiples of 2hkL. 

The atom-optics realization of the kicked rotor con- 
sists in placing laser-cooled atoms in a periodically-pulsed 
standing wave. The pulse (kick) duration r is very short 
compared to the atom dynamics, so that the pulses can 
be assimilated to a Dirac delta function. The Hamil- 
tonian for the center of mass motion of atoms of mass 
M and momentum p along the laser beam propagation 
direction x is thus 
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where we introduced normalized units [1| where time is 
measured in units of the pulse period T, X = 2kLX, 
P = 2kLTp/M, and K = Ahkln^Tr / AlM . With such 
definitions, [X, P] = ik where k = Ahk\T/M plays the 
role of a normalized Planck constant. Note that Pjk ~ 
p/2hkL. 

For A' > 5 the classical dynamics obtained from 
Hamiltonian ([T]) is an ergodic chaotic diffusion j23| which 
leads to a Gaussian momentum distribution whose broad- 
ening corresponds to an average kinetic energy increas- 
ing linearly with time, Ed = Ddt, with Dd ^ In 
the quantum case, after a characteristic localization time 
ti, quantum interference effects lead to a saturation of 
the kinetic energy E to the constant value El = Pl/2 
[24l . [25j . and the average momentum distribution takes an 
exponential shape ^ exp (— |P|/Pl), hence the name of 
"dynamical localization" (DL) given to this phenomenon. 
The instantaneous diffusion coefficient Dn jt) = dE/dt 
starts with an initial value £'o(0) = [26j and tends to 
zero asymptotically. Following we will simply model 
this behavior by an exponential with a characteristic time 
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DL, which is a quantum interference effect, is highly 
sensitive to decoherence and will serve as a decoherence 
probe in the present work. Decoherence effects in the 
kicked rotor have been studied both theoretically [28| and 
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experimentally 0,3 ■ Here, decoherence is essentially due 
to SE, whose probability per kick is 
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where the approximate value on the right corresponds to 
the limit |Al| ^ F, fi. Typical values of 11 in experiments 
are a fraction of lO"^; DL can thus be observed only for 
a hundred kicks or so (2ot . 

Suppose one prepares a sample of cold atoms in a well 
defined momentum state, say P = 0. Interaction with 
the optical potential populates only states of momentum 
Pn = nk, {n integer) and the corresponding wavefunc- 
tion is of the form = J^n'^ne'^''^" \nk). Stimulated 
processes keep well-defined phases and allow the in- 
terference effects responsible for DL. Spontaneous emis- 
sion introduces random phases and puts the atom into 
a mixture of momentum states; in other words, it "re- 
sets" quantum interference, restoring the initial diffusion 
coefficient Dq. The form of the diffusion coefficient D{t) 
with SE is thus a weighted mean of the evolution without 
SE, given by Do{t),ai,nd of the effect of SE events. As 
has been shown in [4I, SE restores diffusion after a time 
t ~ ts{l + Ts)~^ (where Ts = Tits) with an asymptotic 
diffusion coefficient given by 
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Fig.[IJi shows the typical behavior of the kinetic energy 
for increasing levels of SE obtained by numerical simula- 
tion of the quantum dynamics corresponding to Hamilto- 
nian ([1]). Clearly, DL is destroyed and diffusion restored 
as the SE rate increases. The initial momentum distri- 
bution is a square (for simplicity) of width A/? = 0.04^ 
centered at P = 0. For each kick, a Monte Carlo proce- 
dure is used to decide whether a photon is spontaneously 
emitted; if so, the entire momentum distribution is trans- 
lated of a quantity k cos 9, with 9 picked randomly in the 
interval [0,27r), which produces the mixing of quasimo- 
menta. In order to simplify notations, we will measure 
momentum in units of ft, so that, from now on, P„ = n 
and A = An/k. 

Momentum distributions f{P) = \ip{P)\^, are shown 
in Fig. [T|3. For 11 = 0, the distribution is fitted by an ex- 
ponential (red triangles); decoherence effects are clearly 
seen for 11 7^ 0; for instance, the 11 — 0.02 is fitted by a 
Gaussian (black squares). 

Let us now present our method for a transient sup- 
pression of decoherence effects. If no SE event hap- 
pens, the momentum distribution f{P), at any time 
will be nonzero only for momenta in the range A„ = 
[n - A/2, n -I- A/2]. As A < 1, atoms undergoing a SE 
process will mostly populate other momentum classes. 
If, at the end of the kick series, one measures the mo- 
mentum distribution by probing only the A„ momentum 
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Figure 1: (Color online) Effect of spontaneous emission: (a) 
Kinetic energy E{t) for different values of the spontaneous 
emission rate 11; (b) momentum distributions (log scale) at 
t = 500, the n = curve (red triangles) is well fitted by 
an exponential, whereas the 11 = 0.02 curve (black squares) 
is fitted by a Gaussian. Parameters are K = 10, S = 2.9 
{II ~ 12) and the initial momentum distribution is a square 
of width A = 0.04. 
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Figure 2: (Color online) Average kinetic energy with filtering 
(full circles) and without filtering (empty circles) for D = 
0.005 and at iV = 300 kicks (1.5 times the typical SE time, 
— 200 kicks). The filtering reduces the average energy 
at of 60 %. Referring to the saturation value of the kinetic 
energy in absence of SE, the SE effect has been reduced by a 
factor of about 5.5. 



classes, most of the atoms having performed SE will be 
excluded from the resulting signal, which will be effec- 
tively decoherence- free. Experimentally, preparation and 
measurement of f{P) with a precision A ^ 1 can be 
made e.g. using Raman stimulated spectroscopy 2§-3H. 
However, the probability that an atom having performed 
a SE falls back into a detection range A„ increases with 
the number of fluorescence cycles, so this filtering is a 
transient effect. 

Fig. [2] shows the kinetic energy of filtered atoms, E ~ 
E„n'./2 (fuU circles), with P„ ^ f{P)PdP. 
The comparison with the result in the case where there 
is no filtering (empty circles) shows the efficiency of the 
filtering. 

Let us make the reasonable assumption that an atom 
emitting spontaneously a photon has a probability A of 
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falling in one of the ranges A„ and (1 — A) of falling out- 
side any of these ranges. We then consider three popula- 
tions of atoms: Fq, the population of atoms which have 
never performed (at time t) any SE (the momentum of 
these atoms is thus in the ranges A„); Fa, the population 
of atoms in one of the ranges A„ but having performed at 
least one SE, and Fi_a, the population of atoms outside 
any detection range A„. The rate equations for these 
populations are straightforwardly obtained (noting that 
Fa+FA + Fi_A = 1) as: 
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dt 

dFA 
dt 
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The integration of these equations gives: 
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Note that the population of detected (or filtered) atoms 
Fq + Fa= a -|- (1 — A)e~^* is not constant in time and 
tends, as expected, to A as t cxo. 
The kinetic energy of filtered atoms is 



Fa' 
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where Eq and E^ are the total kinetic energy calcu- 
lated over, resp., the populations Fq and F^- The evolu- 
tion of Eq can be obtained from the equation dEo/dt = 
Do{t)FQ — HEq: the first term is the contribution due to 
the kicks [see Eq. ^] and the second term is the deple- 
tion by SE. Thus 



= Dgt, (l 



which has a maximum at time 
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The evolution of E^ is governed by two contributions: 
i) the contribution of the atoms in the population F^ 
which performed at least one SE event and whose kinetic 
energy is due to all past fluorescence cycles [2^, and ii) 
the contribution of the atoms in the population Fq that 
perform a fluorescence cycle at time t and fall into Fa, 
that is 
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dt 
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which gives, after some algebra: 
ADJ, 



Ea = 
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Figure 3: (Color online) (a) Plot of E{t) obtained from the 
numerical simulation (same curves as in Fig. [TJ and best fits 
obtained from Eq. ((9]) with the fit parameters Dq and ts (full 
lines), (b) The 11 — 0.01 curve of plot (a), indicating the 
characteristic times ti « 51 and t2 « 320, and the reduced 
diffusion (dashed line) with coefficient D,- (cf. text). 
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0.005 
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0.02 




25.5 


26.3 


30.7 


36.6 


ts 


50.7 


48.8 


41.3 


32.5 



Table I; Fit parameters Dq and ts corresponding to the curves 
displayed on Fig. |3] 



From Eqs. pUj) and (|13p one obtains an explicit ex- 
pression for E{t) that we do not give here as it somewhat 
cumbersome; let us just consider the limit A, Ts ^ 1 in 
which it can be written: 



E{t) = DqTsIlAe 



(1 + m) 



(14) 



Fig. [HK shows the numerical results for E{t) for different 
values of II. The full lines are flts of the numerical data 
obtained from Eqs. (fTU)) and (I13D . using Dq and tg as 
adjustable parameters (see notes [26l. l27j). whose values 
are given in Table HI 

The main features of the curves in Fig. [3] can be under- 
stood by simple arguments. The curve in Fig. [5]d shows 
two characteristic times: i) time ti [Eq. pd|) ] at which 
the SE starts depleting the population Fq, ii) time t2 
when the contribution of the "incoherent atoms" Fa ex- 
ceeds the "coherent" contribution due to Fg; t2 is thus 
given by the condition F/\{t2) ~ Fo(t2) which leads to 
t2 ~ — lnAn~^. In the absence of flltering, t2 ^ n~^, 
the flltering thus increases this time by a factor | In A|. 

The energy E{t) exhibits a transient reduced diffusion 
regime that corresponds to the quasi-plateau which is 
clearly seen in Fig. the corresponding diffusion co- 
efficient Dr (indicated by the dashed line) can be esti- 
mated by expanding Eq. (|14p around t = 11^^, which 
gives Dr « 2eDqTsARi 2e(l + Ts)ADoo with Doo given 
by Eq. ^ . Thus the condition for the efficiency of filter- 
ing is A < [2e(l + Ts)]~^(~ 0.1 for current parameters). 

Finally, let us consider the effect of flltering process 
on momentum distribution. Fig. |4] displays momentum 
distributions after = 500 kicks, corresponding to the 
curves in Fig. [5^. One sees that in the case 11 — 0.005 
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Figure 4: (Color online) Momentum distributions at A'^ = 500. 
The n = (red triangles) and 11 — 0.005 (blue circles) curves 
are fitted by exponentials, while the 11 — 0.02 (black squares) 
curve is well fitted by a Gaussian. The curve corresponding 
to n = 0.01 (green diamonds) has an intermediate shape. 

(IIA^ = 2.5) the distribution can still be fitted by an ex- 
ponential, showing that DL has survived up to 500 kicks 
(cf. the corresponding kinetic energy curve in Fig. [1)). 
The n = 0.02 [nN = 10) curve is fitted by a Gaussian 
and is thus affected by decoherence. The case 11 = 0.01 
(IITV = 5) has an intermediate shape between exponen- 
tial and Gaussian. 

In conclusion, we have proposed an original method 
allowing the suppression of decoherence effects in the dy- 
namics of the quantum kicked rotor for a finite but quite 
long time compared to the typical duration of current 
experiments. The method can be applied to state of art 
experiments and shall allow more precise studies of sys- 
tems in which decoherence is the main limitation [19, 2(J|. 
Although spontaneous emission was the only source of 
decoherence considered here, it is worth noting that the 
the method is in principle applicable to any source of de- 
coherence which does not conserve quasimomentum (e.g. 
collisions) . 

The authors are grateful to D. Delande for fruitful dis- 
cussions. 
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